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Abstract 

We consider a Hamiltonian theory of spherically symmetric vacuum Einstein 
gravity under Kruskal-like boundary conditions in variables associated with 
the Einstein-Rosen wormhole throat. The configuration variable in the re- 
duced classical theory is the radius of the throat, in a foliation that is frozen 
at the left hand side infinity but asymptotically Minkowski at the right hand 
side infinity, and such that the proper time at the throat agrees with the right 
hand side Minkowski time. The classical Hamiltonian is numerically equal to 
the Schwarzschild mass. Within a class of Hamiltonian quantizations, we show 
that the spectrum of the Hamiltonian operator is discrete and bounded below, 
and can be made positive definite. The large eigenvalues behave asymptoti- 
cally as V2k, where k is an integer. The resulting area spectrum agrees with 
that proposed by Bekenstein and others. Analogous results hold in the pres- 
ence of a negative cosmological constant and electric charge. The classical 
input that led to the quantum results is discussed. 
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I. INTRODUCTION 



One of the most intriguing problems in black hole thermodynamics is the statistical 
mechanical interpretation of black hole entropy. One surmises that black hole entropy should 
reflect an outside observer's ignorance about the quantum mechanical microstates of the 
hole, but it has proved very difficult to characterize what exactly these quantum mechanical 
microstates might be. There are hints that the relevant degrees of freedom may live on the 
horizon of the hole |l]-ffl, in situations where the horizon can be meaningfully defined. There 
is also evidence that black hole entropy describes the entanglement between the degrees of 
freedom in the interior and in the exterior of the hole |§-§| . Recent results from string theory 
H suggest that black hole entropy can be recovered by counting the quantum microstates 
even in situations where the definition of these states presupposes no black hole geometry. 
For reviews, see Refs. P-pH- 

Even prior to Hawking's prediction of black hole radiation ]T2J, the anticipated connection 
between black holes and thermodynamics led Bekenstein |13| to propose that the horizon 
area of a black hole is quantized in integer multiples of a fundamental scale, presumably of 
the order of the square of the Planck length /pi anc k = VfoG* 



c 



-3 



A = akl PlSLnck , (1.1) 

where k ranges over the positive integers and a is a pure number of order one. This proposal 
has since been revived on various grounds; see Refs. [|T4| plf| , and references therein. Although 
the horizon of a (classical) black hole is a nonlocal object, its total area is completely 



determined by the irreducible mass p2f , and one can therefore alternatively view the rule 
( |1 . 1|) as a proposal for the spectrum of the quantum irreducible mass operator .0 As the 
irreducible mass can classically be read off from the asymptotic falloff of the black hole 
gravitational field, one expects such an operator to be sensibly definable even in a quantum 
theory that only refers to observations made at an asymptotically flat infinity. In particular, 
for a Schwarzschild hole, the irreducible mass coincides with the Schwarzschild mass. 



The implications of the area spectrum (O) for macroscopic physics were recently elab- 



orated on by Bekenstein and Mukhanov [26]. Consider, for concreteness, a Schwarzschild 
hole. The area is given in terms of the Schwarzschild mass M by 

A = 16tt f-^^-V M 2 , (1.2) 

\ ^Planck / 

where mpi anc k = \fhcG~ 1 is the Planck mass. Now, ( [Lip and ( |1 . 2| ) imply that M can only 
take discrete values. When the black hole evaporates, it can thus only make transitions 
between the mass eigenstates corresponding to these discrete values. As a consequence, the 
radiation comes out in multiples of a fundamental frequency, which is of the same order 
as the maximum of Hawking's black-body spectrum, and the corresponding wavelength 
is of the order of the Schwarzschild radius of the hole. This means that the radiation 



1 We thank John Friedman and Pawel Mazur for emphasizing this to us. 
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will differ from the black-body spectrum in a way that is, in principle, macroscopically 
observable. For example, if M is of the order of ten solar masses, or 2 x 10 31 kg, then the 
fundamental frequency is of the order of 0.1 kHz, which is roughly the resolving power of an 
ordinary portable radio receiver. The discussion can also be generalized to accommodate a 
nonvanishing angular momentum p"3|| . 

Arguments presented in favor of the area spectrum ( |1 . 1|) include quantizing the 
angular momentum of a rotating black hole ||13||, information theoretic considerations 



|Tj.|T8|.|T^,|26[ , string theoretic arguments |p!q , P!q| , periodicity of Euclidean or Lorentzian time 
]T5| , PT|p7| , |2D| -pi[| , a treatment of the event horizon as a membrane with certain quantum me- 



chanical properties |22| , p4| , and a Hamiltonian quantization of a dust collapse |25| . Recently, 
a membrane model for the horizon |33j recovered an area spectrum that is finer than ( [LID , 
and a calculation within a loop representation of quantum gravity [3|| recovered an area 
spectrum that effectively reproduces the Planckian spectrum for black hole radiation. The 
purpose of the present paper is to give a derivation of (effectively) the spectrum ( |1 . 1| ) from 
a Hamiltonian quantum theory of spherically symmetric Einstein gravity, with judiciously 
chosen dynamical degrees of freedom. 

By Birkhoff 's theorem [|3"5[j , the local properties of spherically symmetric vacuum Einstein 
geometries are completely characterized by a single parameter, the Schwarzschild mass. In a 
classical Hamiltonian theory of such spacetimes, the true dynamical degrees of freedom are 
thus expected to contain information only about the Schwarzschild mass and the embedding 
of the spacelike hypersurfaces in the spacetime. It was demonstrated in Refs. p6|-p9|] that 
this is indeed the case, under certain types of boundary conditions that specify the (possibly 
asymptotic) embedding of the ends of the spacelike hypersurfaces in the spacetime. The 
variables of the reduced theory then consist of a single canonical pair: the coordinate can be 
taken to be the Schwarzschild mass, and its conjugate momentum carries the information 
about the evolution of the (asymptotic) ends of the spacelike hypersurfaces in the spacetime. 
The theory is thus no longer a field theory, but a theory of finitely many degrees of freedom. 
This means, in particular, that quantization of the theory can be addressed within ordinary, 
finite dimensional quantum mechanics. 

Our Hamiltonian theory of spherically symmetric vacuum spacetimes will be built on 
two major assumptions. First, we shall adopt for the spacetime foliation the boundary 



conditions of Ref. p8 |. This implies that the classical solutions have a positive value of the 
Schwarzschild mass, and that the spacelike hypersurfaces extend on the Kruskal manifold 
from the left hand side spacelike infinity to the right hand side spacelike infinity, crossing 
the horizons in arbitrary ways. We shall, however, specialize to the case where the evolution 
of the hypersurfaces at the left hand side infinity is frozen, and the evolution at the right 
hand side infinity proceeds at unit rate with respect to the right hand side asymptotic 
Minkowski time. This means that, apart from constraints, the Hamiltonian will consist 
of a contribution from the right hand side infinity only, and the value of the Hamiltonian 
is equal to the Schwarzschild mass. The physical reason for this choice is that while our 
theory will remain that of vacuum spacetimes, we expect these conditions to correspond to 
physics accessible to an inertial observer at one spacelike infinity, at rest with respect to 
the hole: the proper time of such an observer is our asymptotic Minkowski time, and the 
Arnowitt-Deser-Misner (ADM) mass observed is the Schwarzschild mass. 

Second, we shall adopt as our reduced dynamical variables a canonical pair that is inti- 
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mately related to the dynamical aspects of the Kruskal manifold. Our configuration variable 
a can be envisaged as the radius of the Einstein- Rosen wormhole throat [|32|1 , in a spacetime 
foliation in which the proper time at the throat increases at the same rate as the asymptotic 
Minkowski time at the right hand side infinity. An example of a foliation satisfying these 



conditions can be constructed by taking the Novikov coordinates [^j an d deforming them 
near the left hand side infinity to conform to our boundary conditions there. 

The resulting classical theory has two properties whose physical interest should be em- 
phasized. First, every classical solution is bounded, in the sense that the variable a starts 
from zero, increases to the maximum value 2M, where M is the Schwarzschild mass, and 
then collapses back to zero. This evolution corresponds to the wormhole throat starting 
from the white hole singularity, expanding to the bifurcation two-sphere, and then collaps- 
ing to the black hole singularity. The spacetime dynamics in these variables is therefore, in a 
certain sense, confined to the interior regions of the Kruskal manifold. This property reflects 
the physics observed by an inertial observer at asymptotic infinity, as such an observer sees 
her exterior region of the Kruskal manifold as static. Second, as the proper time on the 
timelike geodesies that form the throat trajectory increases at the same rate as the asymp- 
totic right hand side Minkowski time, one may regard our foliation as a preferred one, by the 
principle of equivalence, for relating the experiences of an inertial observer at the asymp- 
totic infinity to the experiences of an inertial observer at the throat. Note, however, that 
as the total proper time from the initial singularity to the final singularity along the throat 
trajectory is finite, our choice of the foliation implies that the throat reaches the white and 
black hole singularities at finite values of the asymptotic right hand side Minkowski time. 
As the asymptotic right hand side Minkowski time evolves at unit rate with respect to our 
parameter time, this means that the classical theory is incomplete: the classical solutions 
cannot be extended to arbitrarily large values of the parameter time, neither to the past nor 
to the future. 

We quantize the theory by Hamiltonian methods, treating a as a configuration variable, 
specifying a class of 'reasonable' inner products, and promoting the classical Hamiltonian 
into a self-adjoint Hamiltonian operator. For certain choices of the inner product the Hamil- 
tonian operator turns out to be essentially self-adjoint, whereas in the remaining cases the 
class of self-adjoint extensions is parametrized by U(l) and associated with a boundary 
condition at a = 0. We find that the spectrum of the Hamiltonian is discrete and bounded 
below in all the cases. When the Hamiltonian is essentially self-adjoint, the spectrum is 
strictly positive, and in the remaining cases there always exist self-adjoint extensions for 
which the spectrum is strictly positive. A WKB estimate for the large eigenvalues of the 
Hamiltonian yields, via ( [L.2D , the result that the large area eigenvalues are asymptotically 
given byQ 

A ~ 32?rfc /£ lanck + constant + o(l) , (1.3) 

where k is an integer and o(l) denotes a term that vanishes asymptotically at large A. The 
additive constant depends on the choice of the inner product and, when the Hamiltonian 



2 We shall use the symbol ~ to denote an asymptotic expansion throughout the paper. 
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is not essentially self-adjoint, also on the choice of the self-adjoint extension. With two 
particular choices for the inner product, we can verify (and improve on) the accuracy of this 
WKB result rigorously; with another two particular choices, we can rigorously verify the 
accuracy of the leading order term. We can therefore view our theory as producing, from a 
Hamiltonian quantum theory constructed from first principles, the area spectrum with 
a = 327T. 

We shall argue that the discreteness of the quantum spectrum is related to the classical 
incompleteness of the theory. As the variable a classically reaches the singularity at a = 
within finite parameter time, both in the past and in the future, the classical theory can be 
thought of as particle motion on the positive half-line in a confining potential. Whenever 
the self-adjoint Hamiltonian operator is constructed so that the possible 'quantum poten- 
tial' part does not become significant, general theorems guarantee that the spectrum of the 
Hamiltonian will be discrete • I n physical terms, wave packets following classical tra- 

jectories will be reflected quantum mechanically from the origin, and the quantum dynamics 
will in this sense have a quasiperiodic character. In contrast, if the spacetime foliation were 
chosen so that it would take an infinite amount of parameter time for the variable a to 
reach the singularity, the classical theory could be thought of as particle motion on the full 
real line in a potential that is confining on the right but not on the left. In such poten- 
tials, the spectrum of a self-adjoint Hamiltonian operator generically has a continuous part, 
corresponding physically to the fact that wave packets can travel arbitrarily far to the left 
without being reflected. We shall present a simple example of each of these two types of 
foliation, hand-picked so that the Hamiltonian operator becomes easily tractable. From the 
first example we can reproduce the area spectrum ( |1.1| ) with an arbitrary value of the con- 



stant a; in the second example, the spectrum of the Hamiltonian operator will be continuous 
and consist of the full non-negative half-line. 

In addition to the above results for the vacuum theory, we shall also briefly investigate the 
inclusion of a fixed electric charge and a negative cosmo logical constant. With the analogous 
choices for the boundary conditions, the phase space coordinates, and the Hamiltonian 
quantum theory, we again show that the spectrum of the Hamiltonian operator is discrete 
and bounded below. The distribution of the large eigenvalues could presumably be analyzed 
by a suitable generalization of our vacuum techniques; however, we shall not pursue this 
issue here. 

The rest of the paper is as follows. In section || we derive the reduced Hamiltonian theory 
in our phase space variables, starting from Kuchaf's reduced phase space variables and 



performing the appropriate canonical transformation. The theory is quantized in section 11 1 



with considerable parts of the technical analysis deferred to the three appendices. In section 



rV| we discuss the inclusion of the electric charge and a negative cosmo logical constant. 
Section [V| contains a brief summary and a discussion. 

For the remainder of this paper we shall work in natural units, h = c = G = 1. 

II. CLASSICAL WORMHOLE THROAT THEORY 

In this section we present a classical Hamiltonian theory of the Schwarzschild black hole 
in terms of reduced phase space variables that are associated with a wormhole throat, in a 
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sense to be made more precise below. We first briefly recall, in subsection [II A|, Kuchaf 's 



Hamiltonian reduction of spherically symmetric vacuum geometries ||38|| . In subsection [II B 
we derive the throat theory from Kuchaf 's reduced theory via a suitable canonical transfor- 
mation. 



A. Kuchaf reduction 

We start from the general spherically symmetric ADM line element 

ds 2 = -N 2 dt 2 + A 2 {dr + N r dt) 2 + R 2 dtt 2 , (2.1) 
where dfl 2 is the metric on the unit two-sphere, and N, N r , A, and R are functions of t and 



r only. We adopt the falloff conditions of Ref. [38]. These conditions render the spacetime 
asymptotically flat both at r — > oo and r — > — oo, and they make |r| coincide asymptotically 
with the spacelike radial proper distance coordinate in Minkowski space. Each classical 



solution consists of some portion of the Kruskal manifold |3"2" |, such that the constant t 
hypersurfaces extend from the left hand side spacelike infinity to the right hand side spacelike 
infinity, crossing the horizons in arbitrary ways. In particular, the Schwarzschild mass is 
positive for every classical solution. The falloff conditions also guarantee that the four- 
momentum at the infinities has no spatial component: the black hole is at rest with respect 
to the left and right asymptotic Minkowski frames. 

We fix the asymptotic values of N at r — > ±00 to be prescribed i-dependent quantities, 
denoted by N±(t). The Hamiltonian form of the Einstein action, with appropriate boundary 
terms, reads then 

S = j dt J dr (p A A + P R R -NH- N r H r 

- J dt{N + M + + N„M_) , (2.2) 

where H and H r are respectively the super-Hamiltonian constraint and the radial supermo- 
mentum constraint. The quantities M±(t) are determined by the asymptotic falloff of the 
configuration variables, and on a classical solution they both are equal to the Schwarzschild 
mass. We refer to Ref. [Q for the details. Note that when varying the action Q2.2|), N±(t) 



are considered fixed. 

Through a judiciously chosen canonical transformation, followed by elimination of the 
constraints by Hamiltonian reduction, Kuchaf [38] brings the action ( |2.2|) to the uncon- 
strained Hamiltonian form 

S = J dt[pjh-{N + + N_)m] , (2.3) 

where the independent variables m and p are functions of t only. The degrees of freedom 
constitute thus the single canonical pair (m, p): the theory is no longer that of fields, but that 
of finitely many degrees of freedom. The variables take the values m > and — 00 < p < 00, 
and their equations of motion are 
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m = , (2.4a) 
p = -N+ - N_ . (2.4b) 

For our purposes, it will not be necessary to recall the details of the derivation of the 



action (|2.3| ), but what will be important is the interpretation of the reduced theory in terms 
of the spacetime geometry. On each classical solution, the time-independent value of m is 
simply the value of the Schwarzschild mass. By Birkhoff's theorem, m thus carries all the 
information about the local geometry of the classical solutions. The variable p, on the other 
hand, is equal to the difference of the asymptotic Killing times between the left and right 
infinities on a constant t hypersurface, in the convention where the Killing time at the right 
(left) infinity increases towards the future (past). The two terms in the evolution equation 
( p.4b|) arise respectively from the two infinities. Thus, p contains no information about the 
local geometry, but instead it carries the information about the anchoring of the spacelike 
hypersurfaces at the two infinities. Note that p, as the difference of the asymptotic Killing 
times, is invariant under the global isometries that correspond to translations in the Killing 
time. 



B. Hamiltonian throat theory 

We shall now make two restrictions on the reduced theory (|2.3|) . First, we specialize to 
N + = 1 and iV_ = 0. This means that the parameter time t coincides with the asymptotic 
right hand side Minkowski time, up to an additive constant, whereas at the left hand side 
infinity the hypersurfaces remain frozen at the same value of the asymptotic Minkowski time 
for all t. The action reads then 

S = j dt (pm - m) . (2.5) 

The fact that the Hamiltonian now equals simply m reproduces the familiar identification 
of the Schwarzschild mass as the ADM energy, from the viewpoint of asymptotic Minkowski 
time evolution at one asymptotically flat infinity. Our choice as to which of the two infinities 
has been taken to evolve is of course merely a convention, but the choice of completely 
freezing the evolution at the other infinity arises from the requirement that our theory 
describe physics accessible to observers at just one infinity. We shall return to this issue in 
section [V]. 

Second, for reasons that will become transparent below, we confine the variables by hand 
to the range |p| < 7rm. This means that in each classical solution, the asymptotic right hand 
side Minkowski time only takes values within an interval of length 27rm, centered around 
a value that is diagonally opposite to the non-evolving left end of the hypersurfaces in the 
Kruskal diagram. In terms of the parameter time t, each classical solution is then defined 
only for an interval — 7rm < t — to < 7rm, where t = to is the hypersurface whose two 
asymptotic ends are diagonally opposite. 

Consider now the transformation from the pair (m, p) to the new pair (a,p a ) defined by 
the equations 
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v^mfr- 1 - 1 

= V2ma — a 2 + marcsin (1 — a/m) + ^vrm , (2.6a) 
p a = sgn(p) v2ma — a 2 . (2.6b) 

The ranges of the variables are a > and — oo < p a < oo. The transformation is well-defined, 
one-to-one, and canonical. The new action reads 

5 = J dt ( Pa d - H) , (2.7) 

where the Hamiltonian is given by 

H = ~ + a\ . (2.8) 



2 \ a 

The classical solutions are easily written out in the new variables. The value of if on a 
classical solution is just m, and by writing the canonical momentum p a in ( |2.8j ) in terms of 
a and d, one recovers the equation of motion for a in the form 

d 2 = ^-l. (2.9) 
a 

Hence the configuration variable a starts from zero at t = t — 7rm, reaches the maximum 
value 2m at t = to, and collapses back to zero at t = to + 7rm. 

The interest in the variables (a,p a ) is that they have an appealing geometrical interpreta- 
tion in terms of the dynamics of a wormhole throat in the black hole spacetime. To see this, 
we recall that the derivation of the reduced action (|2.3|) from the original geometrodynamical 



action ( |2.2|) in Ref. |38[ relied on the properties of the spacelike hypersurfaces only through 
their asymptotic behavior, but otherwise left these hypersurfaces completely arbitrary. We 
can thus exercise this freedom and seek an interpretation of the variables (a, p a ) in terms of 
a suitably chosen foliation. 

The crucial observation is now that equation ( |2.9|) is identical to the equation of a radial 
timelike geodesic through the bifurcation two-sphere of a Kruskal manifold of mass m, pro- 
vided one identifies a as the curvature radius of the two-sphere and the dot as the proper time 
derivative. This is easily seen for example from the expression of the interior Schwarzschild 
metric in the Schwarzschild coordinates. With these identifications, equations ( |2.6| ) show 
that — p becomes identified with the proper time elapsed along this geodesic from the bi- 
furcation two-sphere, with positive (negative) values of — p yielding the part of the geodesic 
that is in the black (white) hole interior. Thus, if there exists a foliation consistent with 
the falloff conditions of Ref. [fj8|| , intersecting a timelike geodesic through the bifurcation 
two-sphere so that — p agrees with the proper time on the geodesic in this fashion, then the 
quantity a defined by Q2.6a|) is the two-sphere radius along this geodesic. 



It is easy to see that foliations of this kind do exist. Let us briefly discuss an example 
that is closely related to the Novikov coordinates ||32|| . Recall that the geometric idea behind 



the Novikov coordinates consists of fixing a spacelike hypersurface of constant Killing time 
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through the Kruskal manifold, and releasing from this hypersurface a family of freely falling 
test particles with a vanishing initial three-velocity in the Schwarzschild coordinates. The 
coordinates (r, R*) are then defined so that they follow these test particles: the trajectories 
are the lines of constant R*, and on each trajectory r is equal to the proper time. The initial 
hypersurface is r = 0, with R* > and R* < giving respectively the halves living in the 
right and left exterior regions. Now, to arrive at a foliation satisfying our requirements, we 
first deform the Novikov coordinates near the left hand side infinity to accommodate the 
condition iV_ = 0, and we then redefine R* near the infinities in a r- independent fashion so 
as to conform to the radial falloff assumed in Ref. []3S| (the right hand side infinity will then 
have, in the notation of Ref. ||S5| , a falloff with e = 1). The distinguished geodesic through 
the bifurcation two-sphere is given by R* = 0, and the coordinate r agrees by construction 
both with the proper time along this geodesic and with the asymptotic Minkowski time at 
the right hand side infinity. 

Our interpretation of a gives now a geometrical reason for the restriction |p| < 7rm, which 
we above introduced by hand. As a radial timelike geodesic from the initial singularity to the 
final singularity through the bifurcation two-sphere has the finite total proper time 27rm, 
foliations satisfying our conditions do not cover all of the spacetime. The foliations only 
exist for the duration of 27rm in the asymptotic right hand side Minkowski time.Q 

We summarize. Fix a radial timelike geodesic through the bifurcation two-sphere, and 
choose any foliation, consistent with our falloff conditions, such that the proper time along 
the geodesic and the asymptotic right hand side Minkowski time agree on the constant t 
hypersurfaces. Then, the variable a equals the radius of the two-sphere on the distinguished 
geodesic. In particular, if the foliation is chosen so that on each constant t hypersurface, 
the radius of the two-sphere attains its minimum value on the distinguished geodesic, then 
this geodesic and the ones obtained from it by spherical symmetry form the trajectory of 
the Einstein- Rosen wormhole throat. This is the case for example in foliations obtained by 
deforming the Novikov coordinates near the left hand side infinity in the manner discussed 
above. We shall therefore, with a minor abuse of terminology, refer to a as the radius of 
the wormhole throat, and to the theory given by fl2.7|) and ( |2.8|) as the Hamiltonian throat 
theory. 

It is important to note that the configuration variable a is bounded on each classical 
trajectory, reaching the maximum value 2m as the wormhole throat crosses the bifurcation 
two-sphere. This means, in a certain sense, that the spacetime dynamics in terms of our 
configuration variable a is confined 'inside' the hole. This is physically appealing from the 
viewpoint of an observer at infinity: such an observer sees the exterior region of the spacetime 
as static. 



3 Note that although the Novikov coordinates (t,R*) are globally well-defined on the Kruskal 
manifold, the hypersurfaces of constant r extend from one infinity to the other only for |r| < 7rm. 
For larger values of |r|, these hypersurfaces hit a singularity. 
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III. THROAT QUANTIZATION 



In this section we shall quantize the Hamiltonian throat theory of section ||. We saw 
above that the classical Hamiltonian is numerically equal to the Schwarzschild mass, and 
that this Hamiltonian arises as the energy with respect to the Minkowski time evolution at 
one asymptotically flat infinity. The quantum Hamiltonian operator can therefore be viewed 
as the energy operator with respect to an asymptotic Minkowski frame in which the hole is 
at rest. In particular, the spectrum of the Hamiltonian operator becomes the ADM mass 
spectrum of the hole. Our main aim will be a qualitative analysis of this spectrum. 

We take the states of the quantum theory to be described by functions of the configuration 
variable a. The Hilbert space is H := L 2 (R + ; fida), with the inner product 



(3.1) 



where /i(a) is some smooth positive weight function. To obtain the Hamiltonian operator H, 
we make in the classical Hamiltonian (|2.8|) the substitution p a — > —id /(da) and adopt a 
symmetric ordering with respect to the inner product (|3.1|). The result is 



2 



1 d ( fi d 

/j, da \a da 



(3.2) 



For technical reasons, it will be useful to work with an isomorphic theory in which the 
inner product and the kinetic term of the Hamiltonian take a more conventional form. To 



achieve this, we write a 



x 



2/3 



v , and ijj 



mapped to the theory whose Hilbert space is H := L 



(Xl)X2) 



dx X1X2 



The new Hamiltonian operator is 
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_f_ Ax 2 ^ 
dx 2 9 v 



v x- The theory above is then 
f ; dx), with the inner product 

(3.3) 



(3.4) 



where ' = d/dx. By construction, H is a symmetric operator in H . Note that if we had 
retained dimensions, the 'quantum potential' term \v" jv would be proportional to h 2 . 

To completely specify the quantum theory, we need to make H into a self-adjoint oper- 
ator on H . The possible ways of doing this depend on the quantum potential term ^i/'/u. 
For concreteness, we shall from now on take \x = a s , where s is a real parameter. The qual- 
itative results would, however, be analogous for any sufficiently similar /1 with a power-law 
asymptotic behavior at a — > and a — > 00. 

With /i = a s , Hq takes the form 
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d 2 4x 2/3 r(r - 1) 



dx 2 



9 



Xr 



(3.5) 



where 
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r = (2s - l)/6 for s > 2 , (3.6a) 
r = (7-2s)/6 for s < 2 . (3.6b) 

We could have replaced ( |3.6|) by (say) ( |3.6a|) for all s, with (|3.5|) still holding. However, as 



( P-5p is invariant under r — > 1 — r, it will be sufficient to analyze H Q for r > 1/2; this range 
for r is recovered through the definition ( |3.6|) . 

It is easy to see that H has self-adjoint extensions for any r. In the terminology of 



Ref. jy], infinity is a limit point case, whereas zero is a limit point case for r > 3/2 
and limit circle case otherwise (Ref. |41|| , Theorems X.8 and X.10, and Problem 7). For 
r > 3/2, Hq is therefore essentially self-adjoint. For 1/2 < r < 3/2, on the other hand, 
the self-adjoint extensions of H are characterized by a boundary condition at zero and 
parametrized by U(l). 

We now wish to extract qualitative information about the spectrum of the self-adjoint 
extensions of Hq. 

A first observation is that the essential spectrum of every self-adjoint extension of Hq is 
empty (Ref. ||40|| , Theorems XIII. 7.4, XIII. 7.16, and XIII. 7. 17). This means that the spec- 
trum is discrete: the spectrum consists of eigenvalues corresponding to genuine, normalizable 
eigenstates, and the eigenvalues have disjoint neighborhoods. 

Second, we shall show in appendix ^ that every self-adjoint extension of Hq is bounded 
below: the system has a ground state. For r > 3/2, the ground state energy is always 
positive. For l/2<r<3/2, the ground state energy depends on the self-adjoint extension, 
and the situation is more versatile. On the one hand, there is a certain (open) set among 
the self-adjoint extensions within which the ground state energy is positive. On the other 
hand, there exist extensions whose ground state energy is arbitrarily negative. 

Third, we shall show in appendix [FJ that a WKB analysis yields for the squares of the 
large eigenenergies the asymptotic estimate 

-^wkb ~ 2fc + constant + o(l) , (3.7) 

where k is an integer and o(l) denotes a term that vanishes asymptotically at large E. The 
constant depends on r and, for 1/2 < r < 3/2, also on the self-adjoint extension, in a way 
discussed in appendix [B|. We shall also be able to rigorously verify the accuracy of this 
WKB result in the special case r = 5/6, and the accuracy of its leading order term in the 
special case r = 1. 

These properties of the spectrum have consequences of direct physical interest. At the 
low end of the spectrum, the fact that the Hamiltonian is bounded below indicates stability: 
one cannot extract from the system an infinite amount of energy. At the high end of the 
spectrum, the asymptotic distribution of the large eigenenergies yields for the black hole 
area the eigenvalues ( |1.3| ): this agrees with Bekenstein's area spectrum ( |1.1|) , with a = 327r. 



We shall discuss the physical implications of these results further in section M. 



IV. THROAT THEORY WITH CHARGE AND A NEGATIVE COSMOLOGICAL 

CONSTANT 

In this section we shall outline how the throat theory can be generalized to accommodate 
electric charge and a negative cosmological constant. The classical black hole solutions are 
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in this case given (locally) by the Reissner-Nordstrom-anti-de Sitter metric 



35 



, 2 A 2M Q 2 R 2 \ - „ o2 „ 1 . 

1 h — 1 

with the electromagnetic potential one-form 

A = %dT . (4.1b) 
R 

The parameters M and Q are referred to as the mass and the (electric) charge, and the 
cosmological constant has been written in terms of the positive parameter £ as — 3£~ 2 . For 



the global structure of the spacetime, see Refs. [*4"2",f4"3"|. We shall understand the case of a 
vanishing cosmological constant as the limit £ — > oo, and in this case the above solution 
reduces to the Reissner-Nordstrom solution. 

The Hamiltonian structure of the spherically symmetric Einstein- Maxwell system with 
a cosmological constant was analyzed by a technique related to Ashtekar's variables in Refs. 



4*4], 4l"fl . An analysis via a Kuchar-type canonical transformation and Hamiltonian reduction, 



both with and without a negative cosmological constant, was given in Ref. [{46 1. Although 



the focus of Ref. [|6|| was on thermodynamically motivated boundary conditions, which 
confine the constant t hypersurfaces to one exterior region of the spacetime, the discussion 
therein generalizes without essential difficulty to boundary conditions that allow the con- 
stant t hypersurfaces to extend from a left hand side spacelike infinity to the corresponding 
right hand side spacelike infinity, crossing the event horizons in arbitrary ways but crossing 
no inner horizons. The new technical issues arise mainly from the fact that with a neg- 
ative cosmological constant, the left and right infinities are asymptotically anti-de Sitter 
rather than asymptotically flat. The new physical issues arise mainly in the choice of the 
electromagnetic boundary conditions at the infinities. 

We shall here concentrate on the theory where the electric charge is fixed at the infini- 
ties. In the reduced theory, the charge then becomes an entirely nondynamical, external 
parameter, which we denote by q. On a classical solution, q is equal to Q in ([4.11). 

The action of the reduced theory takes the form Q2.3Q . On the classical solutions, m 
is equal to the mass parameter M of ( {4. la| ) . The canonical conjugate p again equals the 
difference in the asymptotic Killing times between the left and right ends of a constant t 
hypersurface. The range of m is m > m crit , where the critical value m crit (q, £) is positive for 
and vanishes for q = 0: this restriction arises from the requirement that the classical 
solutions have a nondegenerate event horizon [|2,4~~j,46[. The range of p is — oo < p < oo. 



The quantities N± determine the evolution of the ends of the hypersurface in the asymptotic 
Killing time. For £ — > oo, N± are simply the asymptotic values of the lapse, whereas for 
< £ < oo they are related to the lapse by a factor that diverges at the infinities. Note that 
in the special case of q = and < £ < oo, we get a theory of vacuum spacetimes with a 
negative cosmological constant. 

Mimicking section |lTj, we freeze the evolution of the hypersurfaces at the left hand side 
infinity by setting AL = 0, and fix the evolution at the right hand side infinity to proceed at 
unit rate with respect to the Killing time by setting N + = 1. The action is given by fl2.5|). 
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The value of the Hamiltonian on a classical solution is then equal to the mass parameter. 
For £ — > oo this reproduces the identification of the mass as the ADM energy from the 
viewpoint of asymptotic Minkowski time evolution at one infinity, just as in the uncharged 
case in section ||. For < £ < oo, we similarly recover the interpretation of the mass 
parameter as the ADM-type energy from the viewpoint of asymptotic anti-de Sitter Killing 
time evolution at one infinity. 

In analogy with Q2.6p, we introduce the new variables (a,p a ) via the transformation 



IP! / ^ , (4.2a) 

1 v/2mH - 1 - q 2 6- 2 - b 2 £~ 2 



p a = sgn(p) y 2ma — a 2 — q 2 — a 4 £ 2 , (4.2b) 

where a + (m, q, £) is the unique positive zero of the right hand side in ( |4.2b| ) for q = 0, and 
the larger of the two positive zeroes for q ^ 0. On a classical solution, a + is the radius of the 
event horizon. To make this transformation well-defined, we again need to restrict by hand 
the range of p. For q = 0, the upper limit for |p| is obtained from ( 4.2a|) with a = 0. For 



q 7^ 0, the upper limit for |p| is obtained from ( [4.2a ) with a = a_, where a_(m, q, £) is the 



smaller of the two positive zeroes of the right hand side in ( 4.2b|) . On a classical solution, 
a_ is the radius of the inner horizon. 

With this restriction on the range of p, the transformation ( f4.2|) is well-defined, one-to- 
one, and canonical. The new action is given by ( |2.7|) with the Hamiltonian 



H= UA + a + £ + *) , (4.3) 



2 \ a a 

and the value of this Hamiltonian on a classical solution is just the mass. 

As in section ||, the theory has an interpretation in terms of a wormhole throat. If there 
exists a foliation with the appropriate falloff conditions [TO , intersecting a timelike geodesic 
through the event horizon bifurcation two-sphere so that — p coincides with the proper time 
from the bifurcation two-sphere along this geodesic, then the quantity a defined by (|4.2a|) 



is the two-sphere radius on this geodesic. If the foliation is chosen suitably symmetric near 
the specified geodesic, we can think of a as the radius of the wormhole throat. We shall now 
examine the existence of such foliations for the different values of the parameters q and £. 

For q = and < i < oo, the classical solution is the Schwarzschild-anti-de Sitter hole. 
The Penrose diagram differs from that of the Kruskal manifold only in that the asymptot- 
ically flat infinities are replaced by asymptotically anti-de Sitter infinities, represented by 
vertical lines [fj^ , |43|l . Foliations of the desired kind clearly exist: the timelike geodesic starts 
at the initial singularity with a = 0, reaches the maximum value of a at the bifurcation 
two-sphere, and ends at the final singularity with a = 0. The situation is thus qualitatively 
very similar to that with the Schwarzschild hole. 

For q ^ and £ —>■ oo, the classical solution is the Reissner-Nordstrom hole with m > 



|q|. The Penrose diagram can be found in Refs. [32,|7|]. Our Kuchaf-type Hamiltonian 
formulation is valid for the part of the spacetime that consists of one pair of spacelike- 
separated left and right asymptotically flat regions and the connecting region that is bounded 
in the past and future by the Cauchy horizons. The solutions to the equations of motion 
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obtained from the Hamiltonian ( |4.3| ) are periodic oscillations in the interval a_ < a < a+, 
but our derivation of this Hamiltonian is only valid on each solution between two successive 
minima of a. Now, it is clear that foliations of the desired kind exist: the timelike geodesic 
starts with a = a_ at the past Cauchy horizon bifurcation two-sphere, reaches a = a + at the 
event horizon bifurcation two-sphere, and ends with a = a_ at the future Cauchy horizon 
bifurcation two- sphere. 

Finally, for q ^ and < i < oo, the classical solution is the Reissner-Nordstrom- 
anti-de Sitter hole, with m so large that a nondegenerate event horizon exists. The Penrose 
diagram is obtained from that of the Reissner-Nordstrom hole by replacing the asymp- 
totically flat infinities by asymptotically anti-de Sitter infinities |42] , |43[| . Our Kuchaf-type 
Hamiltonian formulation is valid for the part of the spacetime that consists of one pair of 
spacelike-separated left and right asymptotically anti-de Sitter regions and the connecting 
region that is bounded in the past and future by the inner horizons. The inner horizons 
are now not Cauchy horizons, as the asymptotically anti-de Sitter infinities render our part 
of the spacetime not globally hyperbolic. As with the Reissner-Nordstrom hole above, the 



solutions to the equations of motion obtained from the Hamiltonian ( |4.3| ) are periodic os- 
cillations in the interval a_ < a < a + , and our derivation of this Hamiltonian is only valid 
on each solution between two successive minima of a. Foliations of the desired kind now 
exist while the timelike geodesic remains sufficiently close to the event horizon bifurcation 
two-sphere. However, it is seen from the Penrose diagram that as the geodesic progresses 
towards the past and future inner horizon bifurcation two-spheres, there will occur a critical 
value of the proper time after which the constant t hypersurfaces would necessarily need to 
become somewhere timelike. Therefore, the throat interpretation can only be maintained 
in the full domain of validity of the Hamiltonian (f4.3|) by appealing to a foliation where the 
constant t hypersurfaces need not be everywhere spacelike. We shall return to this issue in 
section [V]. 

Quantization of the theory proceeds as in section JTJ. The Hamiltonian Hq ( p.5| ) inherits 
the additional terms 



2x 2/3 2 £ 2 

The theorems cited in section |III| show that the additional terms make no difference for the 
existence and counting of the self-adjoint extensions of Hq, and they also show that the 
essential spectrum of any self-adjoint extension of Hq is again empty. For q = 0, the proof 
of the lower bound for the spectrum given in appendix [A] goes through virtually without 
change. For q ^ 0, the charge term modifies the small x behavior of the wave functions, and 
the analysis of the self-adjointness boundary condition is more involved; in particular, there 
is a qualitative difference between the cases 7/6 < r < 3/2, r = 7/6, and 1/2 < r < 7/6, 
arising from whether the next-to- leading term in the counterpart of ve(x) in ( TA^ ) dominates 
the leading order term in the counterpart of ue{x) at small x. However, the modified Bessel 
function asymptotic behavior ( |A6] ) and (XD) still holds, as can be shown by applying the 



series solution method for ue and the 'second solution' integral formula for ve (see, for 
example, Chapter 8 of Ref. The spectrum of every self-adjoint extension is therefore 

again bounded below, and certain self-adjoint extensions are strictly positive. 

One expects that the asymptotic distribution of the large eigenvalues could be investi- 
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gated via a suitable generalization of the WKB techniques of appendix ||. We shall, however, 
not attempt to carry out such an analysis here. 

V. DISCUSSION 

In this paper we have considered a Hamiltonian theory of spherically symmetric vacuum 
Einstein gravity under Kruskal-like boundary conditions. The foliation was chosen such 
that the evolution of the spacelike hypersurfaces is frozen at the left hand side infinity, but 
proceeds at unit rate with respect to the asymptotic Minkowski time at the right hand side 
infinity. The reduced Hamiltonian theory was written in a set of variables associated with 
the Einstein- Rosen wormhole throat: the configuration variable is the radius of the throat, 
in a foliation in which the proper time at the throat agrees with the asymptotic right hand 
side Minkowski time. The classical Hamiltonian is numerically equal to the Schwarzschild 
mass. 

We quantized the theory by Hamiltonian methods, taking the wave functions to be func- 
tions of the classical configuration variable, and including a general power-law weight factor 
in the inner product. The classical Hamiltonian was promoted into a self-adjoint Hamil- 
tonian operator. We found that the spectrum of the Hamiltonian operator is discrete and 
bounded below for all the choices of the weight factor. In the cases where the Hamilto- 
nian operator is essentially self-adjoint, the spectrum is necessarily positive definite; in the 
remaining cases, self-adjoint extensions with a positive definite spectrum always exist. In 
all the WKB estimate gave for the large eigenvalues the asymptotic behavior V2k, 

where k is an integer, and we were able to rigorously verify the accuracy of this estimate 
for four particular choices of the weight factor. The resulting spectrum for the area of the 
black hole agrees with the spectrum ( |1 . 1| ) proposed by Bekenstein and others, with the di- 
mensionless constant a taking the value 32ir. We also showed that analogous results can be 
obtained in the presence of a fixed electric charge and a negative cosmological constant. 

It is perhaps worth emphasizing that the basic postulates of our quantum theory consisted 
of the choice of a Hilbert space and a self-adjoint Hamiltonian operator on it. We did not 
attempt to define more 'elementary' operators, such as those for 'position' or 'momentum', 
self-adjoint or otherwise, from which the Hamiltonian operator could be constructed. This 
issue might merit further study within some geometric or algebraic framework of quantization 

HHED-B 

Even though our theory is that of pure vacuum, our boundary conditions were chosen 
so as to make the results relevant for physics that is accessible to an inertial observer at 
a spacelike infinity. Our spacelike hypersurfaces have evolution at only one infinity, and 
there they evolve at unit rate with respect to the asymptotic Minkowski time. Our classical 
Hamiltonian is therefore the gravitational Hamiltonian with respect to the proper time of 



4 For example, note that the kinetic term of the Hamiltonian operator H ( |3.2|) can be written as 
i/i _1 / 2 p a a _1 /ip a ^ _1//2 , where p a : ip i— > —ifj,~ 1 / 2 (d/da)[i 1 / 2, ip is a symmetric (but not self-adjoint) 
momentum operator that can be regarded as conjugate to the position operator a : if) \—* atp. We 
thank Thomas Strobl for this observation. 
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an inertial observer at the infinity, at rest with respect to the hole. It is thus reasonable 
to think of the eigenvalues of the Hamiltonian operator as the possible outcomes that an 
asymptotic observer would in principle obtain when measuring the ADM mass of the hole. In 
a given (pure) quantum state, the probability for obtaining a given eigenvalue is determined 
by the component of the state in the respective eigenspace in the standard way. Although 
we are here for concreteness using language adapted to a Copenhagen-type interpretation, 
a translation into interpretations of the many-worlds type could easily be made. 

One can also make a case that our throat variable a depicts in a particularly natural way 
the dynamical aspects of the Kruskal manifold. Classically, the wormhole throat begins life 
at the white hole singularity, expands to maximum radius at the bifurcation two- sphere, and 
collapses to the black hole singularity. The dynamics of a is therefore, in a certain sense, 
confined to the interior regions of the Kruskal manifold, and these are precisely the regions 
that do not admit a timelike Killing vector. From a physical viewpoint, using a variable with 
this property is motivated by the fact that an inertial observer at a spacelike infinity sees 
her exterior region of the Kruskal manifold as static. Further, our foliation made the proper 
time at the throat increase at the same rate as the asymptotic right hand side Minkowski 
time: by the principle of equivalence, one may see this as the preferred condition for relating 
the experiences of an inertial observer at the asymptotic infinity to the experiences of an 
inertial observer at the throat. We recall, in contrast, that the reduced phase space variables 
of Refs. [j36|-|38| reflect more closely the static aspects of the Kruskal manifold. Yet another 



set of variables has been discussed in Refs. [ |5^ , |53|1 . 

With this physical picture, the properties we obtained for the spectrum of the Hamil- 
tonian operator acquire consequences of direct physical interest. At the low end of the 
spectrum, the fact that the Hamiltonian is bounded below indicates stability: one cannot 
extract from the system an infinite amount of energy. At the high end of the spectrum, 
in the semiclassical regime of the theory, the discreteness of the spectrum in accordance 
with the area quantization rule ( |1 . 1[ ) yields the macroscopically observable consequences 
discussed by Bekenstein and Mukhanov ||26|| . It should be emphasized, however, that these 
arguments operate at a somewhat formal level, as our theory does not describe how the 
quantum black hole would interact with other degrees of freedom, such as departures from 
spherical symmetry or matter fields. 

On the grounds of classical positive energy theorems, one may feel inclined to exclude by 
fiat quantum theories in which the ground state energy of an isolated self-gravitating system 
is negative. Among our theories, this would amount to a restriction on the self-adjoint ex- 
tension in the cases where the Hamiltonian operator is not essentially self-adjoint. However, 
given the freedom that we have already allowed in the choice of the the inner product, it 
would be a relatively minor further generalization to add to our Hamiltonian operator the 
identity operator with some real coefficient, and to take the coefficient as a new parameter 
in the quantum theory^ The classical limit of the quantum theory would still be correct, 
provided the new parameter is understood to be proportional to Planck's constant. When 
the Hamiltonian operator is not essentially self-adjoint, any given self-adjoint extension can 



5 We thank John Friedman for this observation. 
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then be made positive definite by choosing the new parameter sufficiently large. Note, how- 
ever, that with any fixed value of the new parameter, there still exist self-adjoint extensions 
whose ground state energy is arbitrarily negative. 

In section ||, we obtained the classical Hamiltonian throat theory by first going from 
the geometrodynamical Hamiltonian variables to Kuchar's reduced Hamiltonian theory, and 
then performing a suitable canonical transformation. The interpretation of our variable a 
as the radius of the wormhole throat was only introduced after the fact, by appealing to 
a particular choice of the spacetime foliation. We took the same route in the presence of 
charge and a negative cosmological constant in section [TV]. We have not discussed here how to 
derive a throat theory directly from the unreduced geometrodynamical Hamiltonian theory 
by introducing a gauge and performing the Hamiltonian reduction, but with the appropriate 
gauge choice, the resulting theory should by construction be at least locally identical to ours. 
The only case where one anticipates a difference in the global properties is in the presence 
of both a nonvanishing charge and a negative cosmological constant. In this case, we saw 
in section |TV| that while our throat theory is valid on each classical solution between two 
successive minima of a, the wormhole throat interpretation could be maintained for all 
of this interval only by appealing to a foliation that is not everywhere spacelike. A direct 
Hamiltonian reduction of the geometrodynamical theory in the corresponding wormhole-type 
gauge would thus only yield our theory in a more limited domain, valid on each classical 
solution in a certain interval around a maximum of a. 

Our choice of freezing the evolution of the spacelike hypersurfaces at the left hand side 
infinity was motivated by the desire to have a theory that would describe physics accessible 
to observers at just one infinity. For a vanishing charge and cosmological constant, this 
motivation can be implemented at the very beginning by setting up the Hamiltonian theory 
not on the Kruskal manifold, but instead on the 1LP 3 geon [pl |. To see this, recall 



that the 1LP 3 geon is the quotient of the Kruskal manifold under a freely acting involutive 
isometry: this isometry consists of a reflection of the Kruskal diagram about the vertical 
timelike line through the bifurcation point, followed by the antipodal map on the two- 
sphere. The RP 3 geon has thus only one exterior region, identical to one of the Kruskal 
exterior regions. Further, the MP 3 geon possesses a distinguished MP 2 of timelike geodesies 



through the image of the Kruskal bifurcation two-sphere: in the Penrose diagram 54 



these geodesies go straight up along the 'boundary' of the diagram. The existence of the 
distinguished geodesies reflects the fact that translations in the Killing time on the Kruskal 
manifold do not descend into globally defined isometries of the MP 3 geon. Now, Kuchaf 's 
canonical transformation and Hamiltonian reduction generalize readily to the MP 3 geon |39 



The reduced action is obtained from fl2.3p by setting iV_ = 0, and the momentum p is now 
equal to the difference of the Killing times between the distinguished timelike geodesies and 
the single spacelike infinity. Setting N+ = 1, we are led to the action (|2.5| ). The variable a 
defined by ( |2.6| ) is now equal to the curvature radius of the distinguished MP 2 of timelike 
geodesies. 

The above interpretation of our theory in terms of the MP 3 geon generalizes immediately 
to accommodate a negative cosmological constant. For a nonvanishing charge, on the other 
hand, there exists again an analogous involutive isometry that can be used to quotient 
the manifold, but the electric field is invariant under this isometry only up to its sign. 
This reflects the fact that Gauss's theorem prohibits a regular spacelike hypersurface with 
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just one asymptotic infinity from carrying a nonzero charge. The RP 3 geon interpretation 
does therefore not extend to the charged case with a conventional implementation of the 
electromagnetic field. 

As we have seen, the central input in our classical theory was to parametrize the geometry 
in terms of the radius of the wormhole throat in a judiciously chosen foliation: our variable 
a is the two-sphere radius on a radial geodesic through the event horizon bifurcation two- 
sphere, in a foliation such that the proper time along the distinguished geodesic agrees with 
the asymptotic Killing time at the right hand side infinity. One possible generalization would 
be to relax the requirement that the variable 'live' in the interior regions of the manifold, 
and use instead timelike geodesies that do not pass through the bifurcation two-sphere. To 
examine this, let us for concreteness set the charge and the cosmological constant to zero, 
and let us generalize the canonical transformation (|2.6|) to 



2m[l+(i?S) 2 ] 

db 

= , (5.1a) 
2111b- 1 - [1 + (R* ) 2 ] 



p a = sgn(p) J2ma - a 2 [l + (R*) 2 ] , (5.1b) 



where R$ is a real-valued parameter. The transformation ( |2.6|) is recovered as the special 
case Rq = 0. The ranges of the new variables are again a > and — oo < p a < oo, but 

the restriction for p is now |p| < ?r[l + (-Ro) 2 ] m - The action is given by (|2.7|) with the 
Hamiltonian 



2 



Pl + a 



a ' l + (R* f 



(5.2) 



On a classical solution, the variable a is now equal to the two-sphere radius on a radial 
timelike geodesic whose trajectory is given by R* = Rq, where R* is the Novikov space 



coordinate in the notation of Ref. |32j| . Foliations that make this interpretation possible 
clearly exist. Simple examples are obtained by deforming the Novikov foliation near the left 
hand side infinity as in section jnj, to accommodate the boundary condition AL = 0, and 
(for Rq ^ 0) also near the throat, to prohibit the constant t hypersurfaces from reaching the 
singularity before the geodesic at R* = Rq- Upon quantization along the lines of section 0, 
we find that the spectrum depends on the parameter Rq only through an overall factor: if 

the eigenvalues are denoted by E k , where k ranges over the nonnegative integers, we have 

4^ = [I + TOT 3/4 4 0) • (5-3) 

A much wider generalization would be to relax the requirement, which we above moti- 
vated by the equivalence principle, that the proper time along the throat trajectory agree 
with the asymptotic Killing time. If one allows this freedom, it is not difficult to come 
up with examples of foliations in which the Hamiltonian takes a mathematically simple 
form. As an illustration, let us exhibit two examples in the case of vanishing charge and 
cosmological constant. 
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As the first example, suppose that p is restricted by hand to have the range |p| < 7 1 7rm, 
where 7 is a positive constant. We perform the canonical transformation 

£= v^mcosf^-) , (5.4a) 
7P 



= V 2 /7msin^— J , (5.4b) 

where the ranges of the new canonical variables are £ > and — 00 < p^ < 00. The 
Hamiltonian takes the form 



H = Jh(pl + e) • (5.5) 



We can identify y27£ as the two-sphere radius on a radial geodesic through the bifurcation 
two-sphere in a foliation where the proper time r along this geodesic is 

2m 

f db 

T = - Sgn(P> J V2mt- - 1 ' (5 ' 6) 



27? 

with £ given by ( |5.4a| ). 

To quantize this theory, we adopt the inner product (ipi,ip2) = d£ ipiip2- We define 
the Hamiltonian operator H by spectral analysis as the positive square root of some positive 
definite self-adjoint extension of 'jHsuo, where Hsno '■= \ [~ {d / d£) 2 + £ 2 ~\ is the simple 
harmonic oscillator Hamiltonian operator on the positive half-line. The following statements 
about Hsko can now be verified: (i) the self-adjoint extensions are specified by the boundary 
condition cos(6)ip — sm(8)dil ! /d^ = at the origin, with the parameter 9 satisfying < 9 < 71; 
(ii) the spectrum of each self-adjoint extension is purely discrete; (iii) the eigenfunctions are 
parabolic cylinder functions p5] , and for 9 = and 9 = n/2 they reduce respectively to the 



odd and even ordinary harmonic oscillator wave functions; (iv) if denotes the eigenvalues, 
with k ranging over the nonnegative integers, we have for 9 = and 9 = tt/2 the respective 
exact results = 2k + | and = 2k + and for other values of 9 the asymptotic large k 
expansion ~ 2k + | + ix~ l cot(6>)A;~ 1//2 + o(fc~ 1//2 ); (v) the absence of negative eigenvalues 
is equivalent to the condition that 9 not lie in the interval — 2~ 3 / 2 7r~ 1 [r(l/4)] 2 < tan(#) < 0. 
The resulting spectrum for H therefore agrees asymptotically with the area quantization 
rule ( |1.1| ) for any 9 that makes -ffsHO positive definite. The numerical constant a takes the 
value 327T7. 

As the second example, suppose that p retains the full range —00 < p < 00, and perform 
the canonical transformation 

m 2 

77 = Tp \ ' ( 5 - 7a ) 

cosh I ) 

V2m/ 

^ = sinh(^-) , (5.7b) 

where the ranges of the new canonical variables are rj > and — 00 < p, n < 00. The 
Hamiltonian takes the form 
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(5.8) 



We can identify 2y/rj as the two-sphere radius on a radial geodesic through the bifurcation 
two-sphere in a foliation where the proper time r along this geodesic is 

2m 

f db 

t = -sgn(p) / , (5.9) 

J v2mo _i — 1 

with 7] given by (|5.7a ). 

In the quantum theory we now adopt the inner product (ipijipz) = ' rj~ l drjipiip2- The 
operator H exp := — [r](d/dT])] 2 + rj 2 is essentially self-adjoint, its discrete spectrum is empty, 
and its essential spectrum consists of the non-negative half-line.^ We can therefore define 

1/4 

the Hamiltonian operator H by spectral analysis as (H cxp ) . It follows that the spectrum 
of H is now continuous and consists of the non-negative half-line. 

These examples suggest that the continuity versus discreteness of the spectrum is related 
to the question of whether the wormhole throat reaches the initial and final singularities 
within finite parameter time. There are general grounds to expect this to be the case. 
A classical theory in which the two-sphere radius reaches zero within finite parameter time 
is singular, in the sense that the classical solutions cannot be continued arbitrarily far 
into the past and future. If one quantizes such a theory so that the classical Hamiltonian 
is promoted into a time- independent, self-adjoint Hamiltonian operator, then the unitary 
evolution generated by the Hamiltonian operator remains well-defined for arbitrarily large 
times. If one starts with an initial wave function that is a wave packet following some classical 
trajectory, the quantum time evolution will force the wave packet to be reflected from the 
classical singularity. The reflection is an entirely quantum mechanical phenomenon, and the 
quantum dynamics acquires in this sense a quasiperiodic character. On the other hand, if the 
classical solutions require an infinite amount of time to reach the singularity, one generically 
expects that in the quantum theory a wave packet initially following a classical trajectory 
will just keep following this trajectory, with some spreading, for arbitrarily large times. 
It is clear that these arguments apply without change also in the presence of a negative 
cosmological constant. In the presence of a nonvanishing charge, an analogous discussion 
applies with the singularity replaced by the inner horizon. 

One may hold mixed feelings about a wormhole throat quantum theory that introduces 
a quantum mechanical bounce at a classical singularity or at an inner horizon. On the 
one hand, singularities and inner horizons are places where the classical theory behaves 
poorly, and one anticipates quantum effects to be important. On the other hand, an outright 
bounce may appear an uncomfortably orderly quantum prediction, given that (semi) classical 
intuition associates singularities and inner horizons with collapses and instabilities. Related 



6 These statements follow in a straightforward manner from Chapter VIII of Ref. |4^] after bringing 
H cxp and the inner product to a standard form by the substitution n = e^. Note that if we had 
retained dimensions, this substitution would need to include a dimensionful constant. 
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discussion, in this and related contexts, can be found in Refs. |56"H62]j . While we view 
the model of the present paper as a useful arena where these issues can be addressed in 
relatively explicit terms, the model is undoubtedly dynamically too poor to support confident 
conclusions about the physical reasonableness of a discrete versus continuous black hole 
spectrum. It would be substantially more interesting if our techniques could be generalized 
to models containing degrees of freedom that carry Hawking radiation. 
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APPENDIX A: SEMI-BOUNDEDNESS OF Hi 



In this appendix we shall show that every self-adjoint extension of the Hamiltonian H 
( pj.5| ) is bounded below, and that certain self-adjoint extensions are strictly positive. We 
shall discuss separately the cases r > 3/2, 1/2 < r < 3/2, and r = 1/2. 



1. r > 3/2 

For r > 3/2, H is essentially self-adjoint. Let x be an eigenfunction with energy E. 
With a suitable choice of the overall numerical factor, x is real- valued and has the small x 
expansion 

X(x) = x r (1 + 0{x 2 )) . (Al) 

Both x an d x' are therefore positive for sufficiently small x. If now E < 0, the eigenvalue 
equation shows then that x is increasing for all x > 0. This implies that x cannot be nor- 
malizable, which contradicts the assumption that x is an eigenfunction. Hence the spectrum 
is strictly positive. 



2. 1/2 < r < 3/2 

For 1/2 < r < 3/2, the self-adjoint extensions of Hq form a family characterized by a 
boundary condition at x — 0, and parametrized by U{\). To find these extensions, we note 
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that for any E, the differential equation Hq\ — E% has two linearly independent solutions, 
denoted by ue{x) and ve(x), with the asymptotic small x behavior 



u E (x) = x r (l + 0(x 2 )) , 
Ve ( x ) = x 1 ~ r (1 + 0(x 2 )) . 



(A2a) 
(A2b) 



For real E, both u and v are real- valued. Using the techniques of Ref. [|T|, it is easily 
shown that the eigenfunctions of a given self-adjoint extension of Hq take, up to overall 
normalization, the form 



Xe = cos(9)u E + sm(9)v E , 



(A3) 



where 9 G [0, 7r) is the parameter specifying the self-adjoint extension. The condition (|A3|) 
can be written without explicit reference to the solutions (|A2|) as 







lim 



(2r - 1) cos(^)x r - 1 x - sm(9)x 2{1 ~ r) d ^ X [ 



dx 



(A4) 



We now proceed to obtain a lower bound for the eigenenergies. 

Consider first an extension in the range < 9 < n/2. Comparing the eigenvalue differen- 
tial equation to the corresponding equation with the term |x 2//3 omitted from Hq ( |3.5|) and 
the energy set to zero, one sees that the prospective eigenfunctions with E < are bounded 
below by the function cos(9)x r + sin(^)x 1 ~ r , which does not vanish exponentially at large x. 
However, as the potential increases without bound as x goes to infinity, every eigenfunction 
must vanish exponentially at large x. Hence the spectrum is strictly positive. 

Consider then an extension in the remaining range 7r/2 < 9 < n. Let \ be an eigenfunc- 

1 /2 

tion with energy E < 0. Writing y = (—8E/9) x, the eigenfunction equation reads 



d 2 



dy- 



+ 



r(r — 1) 

y 2 



+ 1 + 



3y_\ 

8E 2 J 



2/3" 



X 



(A5) 



The last term in ( |A5[ ) is asymptotically small at large negative E, uniformly in the interval 
y G (0, M], where M is an arbitrary positive constant. Omitting this last term gives an 
equation whose linearly independent solutions are y 1 / 2 I r _(i/ 2 )(y) and y 1//2 J(i/ 2 )-r(2/), where 
J is a modified Bessel function 
behavior 



of. Therefore x has at large negative E the asymptotic 



X 



(y/2) 1/2 cos(0)r(r+|)(-2£/9) 



-r/2- 



'r-(l/2) 



(y) 



+ sin(^)r(|-r)(-2 J E;/9) ( '- 1)/2 / {1/2) _ r ( 2/ ) 



uniformly for y G (0, M\. The coefficients of the two Bessel functions in 
fixed by comparing the small y expansions of 



(A6) 



have been 
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we can now 



and QA6|) 

By the asymptotic behavior of the Bessel functions at large argument |65| 
choose M so that y 1 ^ 2 I^/2)- r (y) is positive and increasing for y > M/2. For future use, we 
make this choice so that M > 1. Then, the second term in (|A6|) dominates the first term 
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at large negative E, uniformly for M/2 < y < M. Therefore there exists a constant E < 0, 
dependent on r and 9, such that x an d dx/dy are positive at y = M whenever E < E. As 
|r(r — 1)| < 1, equation ( |A5| ) then shows, by virtue of the choice M > 1, that x diverges 
at large y whenever E < E. As x is by assumption normalizable, we thus see that the 
eigenenergies are bounded below by E. 

Note that the lower bound for the eigenenergies is not uniform in 9. For fixed r and 
any given E, there exists a unique self-adjoint extension of Hq such that E is in the spec- 
trum. This is because the differential equation Hqx = Ex has for any E a normalizable 
solution that is unique up to a multiplicative constant, and matching the small x behavior 
of this solution to (|A3|) uniquely specifies the value of 9. One can thus find extensions with 
arbitrarily negative ground state energy. 



3. r = 1/2 

For r = 1/2, the self-adjoint extensions of Hq form again a family parametrized by U(l). 
The boundary condition characterizing the extensions takes the form ( |A3| ), where 9 e [0, it), 
but now with 

u E (x) = x 1 ' 2 (1 + 0(x 2 )) , (A7a) 
v E (x) =u E (x)\nx +0(x 5/2 ) . (A7b) 

Condition ( |A4| ) is replaced by 

= lim | [cos(0) + sin(0) In x] ^ J x ^ - sin(#)aT 1/2 x J . ( A8) 

For the extension with 9 = 0, one sees as above that the spectrum is strictly positive. 
Consider then an extension in the range < 9 < n. Let x be an eigenfunction with 
energy E < 0, and proceed as above. Equation ( |A6| ) is replaced by 

X ~(~!r) / ^ 1/ MH(^)- sin W(7 + |ln(-2 J E/9))]/o(l/)-sin(0)^ o (y)} , (A9) 

where K is the second modified Bessel function |65| and 7 is Euler's constant. At large 
negative E, the term proportional to y 1 ^ 2 Io(y) dominates, and one can argue as above. 
Hence the spectrum is bounded below. 

As in the case 1/2 < r < 3/2, for any given energy E there exists a self-adjoint extension 
such that E is in the spectrum. One can thus again find extensions with arbitrarily negative 
ground state energy. 

APPENDIX B: LARGE EIGENVALUES OF H 

In this appendix we analyze the asymptotic distribution of the large eigenvalues of the 
self-adjoint extensions of H ( |3.5| ). The idea is to match a Bessel function approximation at 
small argument to a WKB approximation in the region of rapid oscillations. We shall again 
discuss separately the cases r > 3/2, 1/2 < r < 3/2, and r = 1/2. 
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1. r > 3/2 



We begin with the case r > 3/2, where Hq is essentially self-adjoint. We shall throughout 
denote by y an eigenfunction with energy E > 0. 

1/2 

Consider first x at small argument. Setting z = (8E/9) ' x, the eigenfunction equation 
reads 



d 2 r(r 
+ 



r 



dz 2 



-1 + 



3z 
8E 2 ) 



2/3' 



X 



(Bl) 



The last term in ( |B"T| ) is asymptotically small at large E, uniformly in the interval z E 
{0,ME^ 2 }, where M is an arbitrary positive constant. Omitting this last term gives an 



equation whose linearly independent solutions are z 1 l 2 J 1 



and z^ 2 N, 



r-(l/2) 



[z), where 



J and N are the Bessel functions of the first and second kind, and only the former solution 



is normalizable at small x ||65|1 . The asymptotic large E behavior of x is therefore 



X oc x 



(1/2) 



(8E) 



1/2 



X 



(B2) 



valid uniformly in any bounded region in x. Here, and from now on, we use the symbol oc to 
denote the asymptotic form at large E, up to a possibly ^-dependent coefficient. Introducing 
two constants Si and 82 that satisfy < 5\ < 62, and using the asymptotic large argument 
behavior of J [B5|, we can rewrite ( |B2p as 



X oc cos 



(8E) 



1/2, 



7ir 

Y 



(B3) 



valid uniformly for 5i < x < 62- 

Consider then the region of rapid oscillations. We take E so large that the eigenvalue 
equation has two turning points, and we denote the larger turning point by Xq. The region 
of rapid oscillations at large E is sufficiently far left of Xq but sufficiently far right of the 
origin. As x decays exponentially right of Xq, the WKB approximation to the wave function 
in the region of rapid oscillations is (see, for example, Ref. |66| ) 



XWKB = \P2(x)] 1/2 COS 



■ rxo 7r 

J dx p 2 {x) - - 



(B4) 



where 



P2 = 3 



2E - x 2/3 



9r(r 



Ax 2 



1/2 



(B5) 



The evaluation of the integral in ( |B4| ) is discussed in appendix |CJ We find 

'(8E) 1/2 x ttE 2 tt 



XWKB OC COS 



o + I + 0(£- 1/2 ) 



(B6) 
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valid uniformly for 5\ < x < 82- Comparing ([B3|j and ( |B6| ) yields for the large eigenenergies 
the WKB estimate 

£wKB~2A; + r + i + o(l) , (B7) 

where k is a large integer, and o(l) indicates a term that goes to zero at large E. 

Note that we have not attempted to control how far the integer k is from the number of 
the eigenvalue that equation (|B7|) is meant to approximate. To obtain a formula that gives 
an approximation to the kth eigenenergy in the limit of large k, one may need to add to the 
right hand side of ([B7|) some even integer. 

As the potential term in the Hamiltonian is smooth without oscillations or step-like 
behavior, and as |£>2P2 2 | vanishes for large E uniformly in 81 < x < 82, one expects the 
WKB estimate to the large eigenenergies to be an accurate one. We shall not attempt to 
investigate the accuracy rigorously for general r, but we shall see below that the accuracy 
can be verified by independent means in the special case r = 5/6, and, to leading order, also 
in the special case r = 1. 

2. 1/2 < r < 3/2 

For 1/2 < r < 3/2, we recall from appendix [A] that the self-adjoint extensions of Hq 
are specified by the boundary condition ( |A3|) at small x. Fixing 9 and proceeding as above, 
we see from the small argument behavior of the Bessel functions [^] that equation ( |B2|) is 
replaced by 

X~(z/2) 1/2 [cos(0)r(r + I)(2E/9)- r/2 J r _ (1/2) (z) 

+ sin(0)r(§ - r)(2E/9f- 1)/2 J {1/2) -r{z) 1 . (B8) 



When 6 = 0, the second term in ( |B8|) vanishes, and we can proceed as above, with xq now 



denoting the larger one of the two turning points for r > 1 and the unique turning point for 



r < 1. The WKB estimate for the large eigenenergies is again given by (|B7|) . When 9 ^ 0, 
on the other hand, the second term in ( |Blj| ) dominates the first term at large E, and r in 
(|B3|) is replaced by 1 — r. The WKB estimate (|B7|) is therefore replaced by 



£wKB~2fc-r-| + o(l) . (B9) 

In the special case r = 5/6, we can verify the accuracy of these WKB results rigorously. 
The eigenfunctions are now x = x 

1/6 U(-±E 2 ,y/2(x 2/3 - E)), where U is the parabolic 



cylinder function that vanishes at large values of its second argument |]55| . The boundary 
condition flA4|) reads 



= cos(6)u(-±E 2 ,-V2E^ - V2sm(6)U'(-lE 2 ,-V2E^j , (BIO) 

where the prime denotes the derivative of U with respect to its second argument. Using 
Olver's asymptotic expansions of parabolic cylinder functions (Ref. |J7|], formula (9.7), and 
the discussion of the derivative on p. 155), we find 
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E 2 ~2k- 2/3 + (E- 8/3 ) for = , (Blla) 
E 2 ~ 2A; + 2/3 + (E^) for = tt/2 , (Bllb) 
£ 2 ~ 2£; + 2/3 + 0(£" 1/3 ) for ^ ^ tt/2 . (Bile) 

This corroborates the WKB results ( |B7"1) and ( p9|) for r = 5/6, and gives an improved bound 
for the error term. 

In the special case r = 1, the theorem in §7 of Ref. |68| yields the rigorous asymptotic 
estimate 

E 2 ~ 2A; + 0(ln£) . (B12) 



(The leading order term of (|B12j ) also follows from Ref. pO |, p. 1614.) This corroborates 
the leading order term in our WKB results ( W^tl arid ( [B9l) for r — I. 

Finally, we note in passing that for r = 5/6, the boundary condition ( |B 1 0| ) and Olver's 
expansions |67j yield the asymptotic relation 



tan(0) ~ - 3l/6 F(l/3)] Q ( { _ Eor ^) (B13) 
for the parameter 9 and the ground state energy E , valid in the limit of large negative E . 



3. r = 1/2 

For r = 1/2, equation ( p8|) is replaced by 

op\ -1/4 

-g-J ^ 1/2 { [oos(e) - sin(^) ( 7 + I In (2E/9))] J (z) + (tt/2) sin(^)iV ( 2 )} , 

(B14) 



where 7 is Euler's constant as before. For any 9, the term proportional to z 1 ! 2 Jq(z) dominates 
at large energies, and equation ( |B3| ) holds with r = 1/2. The WKB estimate for the 



eigenenergies is thus given by (|B7| ) with r = 1/2, for any 9. 



APPENDIX C: EVALUATION OF THE WKB INTEGRAL 



In this appendix we outline the evaluation of the integral in ( [B4| ) at large E. We shall 
throughout assume that E is so large that a classically allowed domain exists (which is a 
restriction only for r > 1), that x is in the classically allowed domain, and that x > 5\, 
where the constant Si was introduced in appendix |B[ 

Returning to the variable a = x 2//3 of the main text, and writing b = (x') 2//3 , the integral 



in the exponent in (B4) takes the form 



"0 



S:= / W2l^Jl-J|=iL , (CI) 
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where ao = x is the (larger) turning point. At large E, a® ~ 2E + 0(E 3 ). 

We fix a constant 5 3 > 0, take E so large that |oo — 2E\ < 83/ E, and restrict a to 
be in the interval 8^ 3 < a < 2E — 8%/E. We can then replace the upper limit of the 
integral in (|C1D by 2E — 83/ E, with the error in S being of order 0(E~ l ). The second 
term under the second square root in ( |(J1[ ) is now uniformly of order 0(E^ 1 ). As the first 
square root is of order 0(E), we can expand the second square root in the Taylor series and 
truncate the series after the third term, with the error in S being of order 0(E~ l ). In the 
truncated integrand, the third term is proportional to b^ u ^ 2 (2E — b)~ 3 ^ 2 , which is uniformly 
of order 0(E^ 3 ^ 2 ), and this term can thus be omitted with the consequence of making an 
error of order 0(E~ l l 2 ) in S. We therefore have 

n2E-S 3 /E 

5~ / db 

■J a 

The integral in flC2|) is elementary. The contribution from the second term turns out to be 
of order 0(E^ 1 ^ 2 ), and in the first term the upper limit can be replaced by 2E with the 
consequence of making an error of order 0(E~ l ). We thus obtain 

S ~ \{E - a)V2Ea - a 2 + \E 2 arcsin (1 - a/E) + \irE 2 + 0(E- 1/2 ) . (C3) 

Finally, we restrict a to be in the interval 8^ 3 < a < 8 2 ^ 3 , where the constant 82 was 
introduced in appendix |B|. A large E expansion of ( |C3|) then gives 

which leads to (|B6D. 



V2Eb -b 2 - 



9r(r — 1) 
8b 2 V2Eb - 6 f 



+ 0(£- 1/2 ) 



(C2) 
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